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OBJECTIVES 


This  study  was  part  of  the  first  phase  of  ONR’s  Mine  Burial  Prediction  program  (Coastal  Geosciences, 
Code  321).  The  original  objectives  of  this  study  were  to: 

•  Develop  a  numerical  procedure  for  modeling  mine  burial  as  a  time-dependent  process  with  a 
focus  on  sub-seabed  liquefaction  and  sediment  mechanics. 

•  Collect  necessary  material  parameters,  including  pore  pressure  generation  characteristics  under 
cyclic  loading  conditions,  from  planned  field  and  laboratory  tests. 

•  Integrate  findings  from  this  study  with  parallel  modeling  efforts  underway  as  part  of  ONR’s 
Mine  Burial  Program  (MBP),  which  is  investigating  burial  processes  due  to  bedform  migration 
and  local  scour. 

Funding  for  the  project  was  terminated  after  the  first  two  years  of  an  anticipated  5-year  program  due  to 
a  drastic  reduction  in  funds  for  the  overall  Mine  Burial  Prediction  program,  requiring  that  some 
projects  be  dropped.  Therefore,  substantial  progress  was  only  achieved  on  the  first  objective  and  the 
laboratory  component  of  the  second  one.  No  major  field  experiments  had  been  completed  by  the  time 
funding  for  this  project  was  ended. 


OUTCOMES 

1.  Research. 

A  methodology  and  an  associate  numerical  code  (MINE-B)  were  developed  for  the  analysis  of 
liquefaction  in  shallow  sediments  due  to  cyclic  loading  by  surface  water  waves.  In  addition,  a 
laboratory  testing  program  was  conducted  on  specimens  collected  from  the  Tampa  Bay 
experiment  site  to  develop  a  more  accurate  non-linear  cyclic  stress-strain  model  for  use  in 
MINE-B  or  other  similar  numerical  codes.  Scientific  achievements  are  described  in  the 
publications  listed  below  and  more  recent  findings  are  highlighted  in  this  report. 

2.  Publications  &  Conference  Presentations. 

Brandes,  H.G.  and  Riggs,  H.R.  (2002a).  Modeling  of  seabed  liquefaction  and  other  processes 
responsible  for  mine  burial.  Proceedings  of  OMAE’02,  21st  International  Conference  on 
Offshore  Mechanics  and  Arctic  Engineering,  OMAE2002-28242,  pp.  1-8. 

Brandes,  H.G.  and  Riggs,  H.R.  (2002b).  Numerical  modeling  of  seabed  liquefaction  due  to 
surface  water  waves.  Proceedings  of  the  12th  International  Offshore  and  Polar 
Engineering  Conference,  2:558-562. 

Brandes,  H.G.,  Seidman,  J.  &  Wang,  S.  (2005).  Cyclic  simple  shear  testing  and  modeling  of 
two  calcareous  sediments.  15th  International  Offshore  and  Polar  Engineering  Conference, 
Seoul,  Korea,  June  19-24. 

3.  Students. 

Shentang  Wang  -  Ph.D.  Programming  and  constitutive  modeling.  Mr.  Wang  is  expected  to 
complete  his  dissertation  work  in  2005.  Title:  OpenSees  model  development  for  static  and 
cyclic  soil-structure  interaction  modeling. 


Jason  Seidman  -  B.S.  Assisted  with  equipment  development  and  laboratory  testing  of  field 
specimens.  Mr.  Seidman  is  expected  to  continue  with  graduate  studies  upon  graduation  in 
2005. 

4.  Equipment. 

A  direct  simple  shear  system  was  acquired  with  funds  from  this  project  and  a  matching 
contribution  from  the  University  of  Hawaii.  This  major  piece  of  instrumentation  enabled  us  to 
carry  out  cyclic  laboratory  tests  to  simulate  loading  from  surface  water  waves.  Results  from 
these  tests  are  presented  by  Brandes  et  al.  (2005)  and  are  briefly  discussed  in  this  report.  The 
equipment  adds  a  unique  state-of-the-art  laboratory  testing  capability  to  the  geotechnical 
laboratory  at  the  University  of  Hawaii.  It  continues  to  be  used  extensively. 

5.  Research  Meetings  &  Presentations. 

ONR  Mine  Burial  Meeting ,  March  5,  2002  -  Phoenix. 

2nd  Annual  ONR  Mine  Burial  Prediction  Workshop ,  January  15-17,  2002  -  San  Diego. 

Coastal  Processes  Driving  Mine  Burial  Workshop ,  January  30-February  1,  2001  -  St. 
Petersburg. 


NUMERICAL  AND  CONSTITUTIVE  MODELING 

A  finite  element  code  based  on  Biot’s  physically  consistent  theory  of  consolidation  for  poro- 
elastic  media  was  developed  in  an  effort  to  predict  time-dependent  mine  burial  due  to  the  effects  of 
surface  water  wave  propagation.  Initially,  a  linear  elastic  constitutive  sediment  model  was 
implemented  for  comparison  of  2D  stress  and  strain  predictions  with  closed  form  solutions  (Yamamoto 
et  al.,  1978;  Madsen,  1978).  A  simple  case  study  was  selected  for  this  purpose.  It  consisted  of  a 
uniform  seabed  of  linearly  elastic  sand  with  a  single  sinusoidal  surface  water  wave  propagating  along 
the  sea  surface.  Linear  wave  theory  was  assumed.  Predictions  obtained  with  MINE-B  for  a  level 
seabed  agree  quite  well  with  theory  (Figures  1-2;  Brandes  &  Riggs,  2002). 

Although  this  problem  is  overly  simplistic,  it  does  allow  for  model  verification  and  can  be  used 
for  preparing  preliminary  liquefaction  charts  for  estimating  the  extent  of  momentary  liquefaction  as  a 
function  of  water  depth  and  wave  characteristics.  The  depth  to  which  liquefaction  is  predicted  may 
correlate  with  mine  burial,  since  by  definition  a  liquefied  seabed  cannot  support  any  object  whose 
weight  is  larger  than  that  of  the  sediment.  An  example  of  such  a  chart  is  shown  in  Figure  3,  which  can 
be  used  to  predict  the  depth  of  momentary  liquefaction  below  the  seabed  that  occurs  under  the  crest  of 
single  linear  waves.  The  conventional  criteria  for  the  onset  of  liquefaction  proposed  by  Seed  and  Lee 
(1966)  is  used  here: 


-fsl.O  (1) 

yz 

In  other  words,  this  criteria  assumes  that  liquefaction  will  occur  when  the  ratio  of  excess  pore  pressure 
p  exceeds  the  effective  overburden  stress  y’z.  Here  we  have  assumed  that  y  -10  kN/m3  and  z  is  the 
distance  below  seabed. 
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Figure  1.  Excess  pore  pressure  in  2-D  seabed  deposit  after  150  wave  cycles  (wave  length=48  m, 
period=6  seconds,  wave  amplitude=2.4  m) 


Free  Field  Response 


Figure  2.  Calculated  versus  predicted  pore  pressure  ratio  below  crest  of  wave. 
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Figure  3.  Predicted  depth  of  liquefaction  below  seabed  for  single,  linear  surface  wave  (d=water  depth, 
L=wave  length,  H=wave  height.  Period  T  is  given  by  linear  wave  theory). 


A  linear  elastic  constitutive  model  is  not  adequate  when  strains  and  stresses  induced  in  the 
seabed  exceed  a  certain  level  and  the  material  begins  behaving  in  a  non-linear  and  irreversible,  i.e 
plastic,  fashion.  In  order  to  account  for  such  effects,  work  was  begun  on  implementing  a  more  robust 
non-linear  cyclic  constitutive  model  (Brandes  and  Wang,  2004).  This  model  takes  into  account 
gradual  stiffness  and  strength  degradation  with  repetitive  cycles  of  loading  and  allows  for  pore 
pressure  dissipation  due  to  pore  fluid  diffusion.  An  example  of  cyclic  shear  stress  and  pore  pressure 
predictions  using  this  model  is  shown  in  Figure  4.  It  is  clearly  seen  that  liquefaction  develops  near  the 
end  of  the  test,  resulting  in  very  large  strains  and  hence  failure.  The  computations  used  a  set  of 
reasonable  but  assumed  material  parameters.  In  order  to  adapt  the  model  for  use  with  sediments  from 
the  field,  a  laboratory  testing  program  was  conducted  on  specimens  from  the  Tampa  Bay  field 
experiment  site.  This  is  discussed  in  the  next  section. 
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LABORATORY  TESTING 

The  laboratory  testing  program  was  designed  to  provide  material  parameters  for  the  constitutive  model. 
It  was  conducted  on  fine  carbonate  sand  collected  during  the  Mine  Burial  program’s  Tampa  Bay  field 
experiment.  A  summary  of  tests  conducted  to  date  is  shown  Table  1  and  further  details  on  results  and 
analysis  are  given  in  Brandes  et  al.  (2005).  An  example  of  the  shear  stress-strain  response  is  illustrated 
in  Figure  5.  Note  the  stiffness  degradation  that  occurs  as  the  test  progresses,  much  as  the  model 
predicts  (Figure  4,  top  graph).  One  important  difference  between  model  and  test  results  is  that 
permanent  shear  strains  accumulate  toward  the  end  of  the  test  as  the  material  approaches  liquefaction. 
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It  is  therefore  evident  that  the  constitutive  model  needs  to  be  modified  to  allow  for  irreversible 
deformations.  On  the  other  hand,  the  model  predictions  in  Figure  4  do  account  for  the  dual  stiffening 
and  softening  behavior  within  a  single  cycle  that  are  apparent  toward  the  end  of  the  test. 


Bender  elements  were  used  during  testing  to  determine  stiffness  at  very  low  strains,  which  complement 
the  stiffness  data  at  larger  strains  from  the  cyclic  testing.  An  example  of  the  results  is  shown  in  Figure 
6  and  more  details  are  presented  in  Brandes  et  al.  (2005). 

Table  1.  Cyclic  simple  shear  test  program 


Test  designation 

Normal  stress 
(kPa) 

Density 

(gm/cm3) 

Cyclic  stress  amplitude 
(kPa) 

SS-TB-LD-1 

86 

1.47 

11.4 

SS-TB-LD-2 

143 

1.46 

11.4 

SS-TB-LD-3 

200 

1.39 

11.4 

SS-TB-LD-4 

257 

1.40 

11.4 

SS-TB-HD-1 

86 

1.76 

11.4 

SS-TB-HD-2 

143 

1.79 

11.4 

SS-TB-HD-3 

200 

1.73 

11.4 

SS-TB-HD-4 

257 

1.66 

11.4 

CL 

jc 


C/3 


Shear  Strain  (%) 

Figure  5.  Cyclic  stress-strain  behavior  for  test  SS-TB-HD-4. 

SUMMARY 

Significant  progress  has  been  made  in  understanding  the  cyclic  stress-strain  behavior  of  granular 
calcareous  sediments  through  laboratory  testing  to  guide  the  development  of  an  appropriate 
constitutive  model.  A  specific  model  has  been  developed  to  the  point  where  it  is  now  ready  for 
implementation  in  a  numerical  code.  Once  this  is  achieved,  burial  of  objects  into  the  seafloor  can  be 
computed  in  the  conventional  way  by  creating  2D  or  3D  mine-seabed  finite  element  geometries  and 
conducting  a  time-dependent  analysis. 
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Figure  6.  Shear  wave  measurement  using  bender  elements.  Test  SS-TB-HD-4. 
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